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Infrared singularities in massless gauge theories are known since the foundation of quantum field 
theories. The root of this problem can be tracked back to the very definition of these long-range 
interacting theories such as QED. It can be shown that singularities are caused by the massless 
degrees of freedom (i.e. the photons in the case of QED). In the Bloch-Nordsieck model the absence 
of the infrared catastrophe can be shown exactly by the complete summation of the radiative cor¬ 
rections. In this paper we will give the idea of the derivation of the Bloch-Nordsieck propagator, 
that describes the infrared structure of the electron propagation, at zero and finite temperatures. 
Some ideas of the possible applications are also mentioned. 


I. INTRODUCTION 

The infrared (IR) limit of quantum electrodynamics (QED) is known to be plagued by singularities caused by the 
photons. This phenomenon is known as the infrared catastrophe, and it can be found in any quantum held theory 
(QFT) which involves massless helds. The development of QFTs started around 1930 with QED, therefore, in most of 
the cases the subjects of the computations were electromagnetic quantities. The methods used for the calculations were 
mostly the direct extension of the PT from quantum mechanics. Physicists back then, who were doing computations 
in QED, immediately faced IR divergences when calculating hrst order perturbative corrections to the Bremsstrahlung 
process, due to the low frequency photon contributions. The core of the problem lays in the fundamental dehnition 
of QED, namely, that we assume the existence of a free theory, i.e. the existence of asymptotic states. However, 
such states are difficult to define in a theory where we have long-range interactions. As a consequence, one cannot 
truly define the asymptotic states described by the Fock representation of free theory Hilbert space, on which the PT 
is performed. Thus, we need to search for a non-perturbative solution to prevent these difficulties. An alternative 
approach to this problem was provided by Bloch and Nordsieck in 1937 in their remarkable work on treating the 
infrared problem [1]. The divergencies are caused by the fact that in a scattering process an infinite amount of long 
wavelength photons are emitted, and these low energy excitations of the photon field are always present around the 
electron in the form of a ’’photon cloud”. This shows us essentially that the observed particle is in fact very different 
from the one we call the bare particle: they can be considered as dressed ’’quasi particle” objects whose interactions 
cannot be described through PT entirely. In this paper, we will show the emergence of the infrared catastrophe and 
then we will introduce the Bloch-Nordsieck (BN) model, which was designed in order to imitate the low energy regime 
of QED. We will discuss the breakdown of the PT due to the IR catastrophe, however, it is possible to obtain the 
exact full solution by using the Ward-Takahashi identities embedded into the Dyson-Schwinger (DS) equation. 


II. THE INFRARED CATASTROPHE 

The easiest way to demonstrate the IR catastrophe is the following. Suppose that E{uj) is a finite amount of 
electromagnetic energy emitted by an accelerated charge in the frequency band [u),u! duj]. Each photon carries an 
energy of hw, hence the average number of emitted photons in this band is h = E{Lo)/huj. If we take the limit w —)■ 0 
the average number of photons will diverge provided that limt,j_i.o E{uj) ^ 0 (which is fulfilled). Thus we can see that 
an infinite number of soft photons are present at any scattering process (c.f. [2-4]). In the following we will apply 
the convention used by the particle physics community, i.e. = c = 1 in the further computations. 

We should get the same result using quantum computations, however, relying on PT gives different result: already in 
the first order of the PT the probability of emitting one soft photon in a scattering process will diverge logarithmically 
[2, 3]. This is in complete contradiction that we have found with the semi-classical line of thought a few lines above. 
In fact, it turns out it is not enough to take into account the tree level diagrams but we will also need to include the 
virtual corrections to cancel out these infrared divergencies. This can be done in all orders of PT and by summing 
up these corrections (combining real and virtual corrections) to infinite order we can obtain a well defined probability 
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measure. More precisely, it can be shown (c.f. [4]) that the probability of emitting n soft photons in the process has 
the form 


PAE,u:^iAn) = l(/|5'b)|^ ^ 

nl 


2a -q^ , /\E 

-ill - 

StT TTl ^min 


-^^In - 


( 1 ) 


Here the first factor is the absolute square of the scattering amplitude without radiation emitted, a is the fine 
structure constant (a = e^/ 47 r, with e being the electric charge); vr? and (f' is the electron mass and the transferred 
four momentum q = p — p' (and —q^ > 0), respectively. There are two energy scales that have been introduced: Wmm 
and AE. The former is an artificial IR regulator (’’mass” for the photon field) and the second is the resolution of the 
detector that performs the measuring in the process: photons having energy lower than AE are not being detected 
at all. Hence, we only consider the interval where the photons energy are uj € [ujmim AE], However, (1) will give 0 
for any finite n while taking the artificial mass of the photon to zero as we should: 


lim (n) = 0 . 

OJmin^O 


( 2 ) 


This means that the probability of emitting any finite number of soft photons during the scattering process is zero. 
On the contrary: if we perform a summation over all possible photon numbers that can be emitted we will get a finite 
result ^ 

°° _ 2 

P{AE) = J2PAE{n) = \{p'\S\p)fe-"-^^'^^. (4) 

n—0 

As we can see, indeed, we obtained a finite probability for the process of infinite emitted photons from the quantum 
computation. However, we still need to keep the sensitivity of the detector finite in order to get this result. Although 
theoretically it is possible to take the limit AE —)■ 0, however, in reality it will never happen since there is no such as 
a detector with perfect resolution. 


III. THE SOFT PHOTON CONTRIBUTION TO THE ELECTRON STRUCTURE: THE 

BLOCH-NORDSIECK MODEL 

The BN model was made to give an insight in the analytic structure of the cancellation of the infinities in the IR 
regime. Investigating this model will lead us to the exact result of the propagator of the fermion which is surrounded 
by a cloud of photons. Although the result is known long ago [1], and it can be considered as a textbook material [5], 
we will sketch the derivation used in [ 6 ] which can be extended to finite temperatures the most easily [7]. The main 
idea of the BN model, that simplifies the theory tremendously, is to replace the gamma matrices 7 ^ by a four-vector 
Uf^ that can be considered as the four-velocity of the fermion, and hence the fermion field is represented by a scalar 
field. This simplification is well justified in the IR regime: the soft photons which take part in the interaction will not 
have enough energy for pair production, moreover, not even enough to flip the spin of the electron. It implies that 
the photon propagator will not have any corrections, i.e. the exact photon propagator is the free one in this model. 
Thus the Lagrangian reads 

C = + -m)^, iD^ = id^- eA^, (5) 

where 4* and A^ is the fermion and photon field, respectively. The fermionic part of the Lagrangian is Lorentz- 
covariant, therefore we can relate the results with different choice by Lorentz transformation. This makes possible 
to work with u = (uq, 0,0,0) without loss of generality. In fact, we can perform a Lorentz-transformation where 
Au = (uq, 0, 0, 0). Since is a four-velocity then uq = 1; if it is of the form u = (1, v), then it is uq = Vl — v^. After 
rescaling the field as 4* —>■ 4*/and the mass as to —>■ ugm, the Lagrangian reads 

£ = -^F^,E^'' + ^(Wo-m)^. (6) 
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This result can be found in [4] where functional techniques are used, however, in [2] the following formula is given for the differential 
cross section: 



dQ \ dQ J Q 

Here the first factor corresponds to the hard scattering and in the exponent we can find the famous Sudakov double logarithm. 
difference between the results in [2] and [4] originates from the different approximations that are used. 


(3) 


The 
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Using this reference frame simplifies the calculations [6]. If necessary, the complete u dependence can be recovered 
easily, however, at finite temperatures we need to use numerics if we want to switch to another reference frame due 
to the lack of the Lorentz symmetry. 


A. The Bloch-Nordsieck model at T = 0 

In order to derive the dressed fermionic propagator, we are going to use a system of self-consistent equations which 
involves the DS equation, Dyson’s series and the Ward-Takahashi identities. We will begin with the DS equation; it 
reads in real and momentum space, respectively: 

T.{x — y) =—ie^ J d'^w j d!^zQ{x — w)u^G^^{x — z)T^{z-,w,y), 

S(p) = -*e^y ^(P - k)u>^G^^ik)T‘'{k;p- k,p), (7) 

where is the photon propagator, Q is the fermion propagator and is the vertex function. The diagrammatic 
representation of equation (7) can be seen in Fig. 1. The DS equation describes the self-energy of the fermion. To ob¬ 
tain the full expression, we need to treat Q as the exact fermion propagator and keep the photon propagator undressed 
(i.e. on tree-level), which coincides with the exact one in the framework of the BN model. The vertex correction is 
composed from both propagators but it can be simplified using the Ward-Takahashi relations, as we will see. 



FIG. 1: The diagrammatic representation of the Dyson-Schwinger equation. The double line is for the dressed 
fermion, the wavy line is for the photon propagator. The black blob denotes the full vertex function. The bold 
letters are for space-time points and the Greek letters denote the Lorentz indices. At finite temperatures we also 
need to consider the regular letters which are the Keldysh indices. 


Now that we have the formula for the self-energy, we need another equation which expresses the fermion propagator 
as a function of the self-energy. For this purpose Dyson’s series can be used [2]. It can be shown that summing up all 
the radiative corrections to the free propagator will lead us to a sum of a geometric series (see Fig. 2). This results 
in an algebraic form: 

Q~^{p) = Z{p-m)-Y,{p), (8) 


where we introduced Z as the wave function renormalization factor coming from the vertex part. At this point 
we are not done yet, since in order to compute anything from (7) and (8) we would need an approximation of the 
vertex function F'^ to close the system of equations. However, in the case of the BN model we can relate the vertex 
correction to the fermion propagator in a unique way. This is the key point of the computation presented in [6], and 
it can be achieved using the Ward-Takahashi identities. This exact equation is coming from the current conservation 
analogously to the QED case [2, 3]: 

kfj,T^{k;p- k,p) = Q~^{p) - Q~^{p- k). (9) 


In this model, however, the vertex function is proportional to . In principle, the Lorentz-index in this model can 
come from or from any of the momenta. But, since the fermion propagator depends on the four-momentum in the 
form u^p^, the fermion-photon vertex does not depend on the momentum components which are orthogonal to u^. 
Therefore the Lorentz-index which comes from k^ in fact comes from the longitudinal part of k^, i.e. proportional to 
So we can write T^^{k]p— k,p) = u^r{k]p— k,p). Using equations (7), (8) and (9) we can compute the exact 
fermion propagator. The detailed calculation can be found in [6], where the Feynman gauge is used and the UV 
renormalization is performed systematically. This all boils down to the to the following differential equation: 



TT 


{up - m)Q' + Q 


( 10 ) 
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and its solution is: 


Q (p) oc {up — m) ^ 


( 11 ) 


^ = —>— + ^ ^ 0 ^ 0 ^^ 
= —>— + — 


FIG. 2: The diagrammatic representation of the Dyson’s series. The figure shows that the exact fermion propagator 
can be obtained as the sum of the radiative corrections. The double line represents the full fermion propagator, the 
grey blob is the self-energy which incorporates the radiative corrections. It can be shown that this sum is a 

geometric series and has the algebraic form of (8). 

This is indeed the exact solution for the fermion propagator in the BN model [1]. Comparing it to the free 
propagator, which has the form of l/{up — m) the dressed solution has a power-law correction. This correction is 
due to the fully quantised soft radiation field, which shows in this respect that the electron cannot be thought as an 
independent object, but it is rather quasiparticle surrounded by the soft photon cloud. In the following we are going 
to give the solution for the finite temperature case. 


B. The Bloch-Nordsieck model at T > 0 

In this section we will show the solution for the BN model at finite temperatures. There are a few studies on the 
finite temperature BN model with different approximations [8-12], however, the only analytically closed formula for 
the spectral function of the fermion is derived in [7]; we will present the results from the latter. We already discussed 
that the parameter as the fixed four-velocity of the fermion implies that the Bloch-Nordsieck model describes the 
regime where the soft photon fields do not have energy even for changing the velocity of the fermion (no fermion 
recoil). This leads to the interpretation that the fermion is a hard probe of the soft photon fields, and as such it is 
not part of the thermal medium [8, 9]. 

We are interested in the finite temperature fermion propagator. To determine it, we use the real time or Keldysh 
formalism (for details, see [13]). Here the time variable runs over a contour containing forward and backward running 
sections (Ci and (72). The propagators are subject to boundary conditions which can be expressed as the KMS (Kubo- 
Martin-Schwinger) relations [13]. The physical time can be expressed through the contour time t = T{t). This makes 
possible to work with fields living on a definite branch of the contour, 4'a(t,x) = 4'(Ta,x) where T{Ta) = t, and 
Ta G Ca for a = 1,2; and similarly for the gauge fields. The propagators are matrices in this notation, in particular 
the fermion and the photon propagator, respectively, reads as: 

iGabix) = (Tc^a{x)'^l{0)), 

lG^,,ab{x) = {TcAf,a{x)A,i,{0)) , ( 12 ) 

where Tc denotes ordering with respect to the contour variable (contour time ordering). For a generic propagator 
Gii corresponds to the Feynman propagator, and, since the G2 contour times are always larger than the Gi contour 
times, G 21 = G^ and G 12 = G^ are the Wightman functions. The KMS relation for a bosonic/fermionic propagator 
reads Gi2{t,x) = zLG2i{t — i/ 3 ,x) which has the following solution in Fourier space (with k = (A:°,k)) 

iGi2{k) = ±n±{ko)g{k), 

iG2i{k) = (1 ± n±){ko)g{k), (13) 

where 


n±{ko) = ^ and g{k) = iG 2 i{k) - iGi 2 {k) (14) 

are the distribution functions (Bose-Einstein (-I-) and Fermi-Dirac (-) statistics), and the spectral function, respectively. 
It is sometimes advantageous to change to the R/A formalism with field assignments 4'i_2 = 4'^ ± 4*0/2. Then, one 
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has Gaa = 0 for both the fermion and the photon propagators. The relation between the propagators in the Keldysh 
formalism and the R/A propagators reads 

Grr = -^, Gll = Gra + Gl2, Q = *Gra ~ iGar- (15) 

The Gra propagator is the retarded, the Gar is the advanced propagator, Grr is usually called the Keldysh propagator 
in the framework of the R/A formalism (not to confuse with the propagators in the Keldysh formalism). 

At zero temperature, as we could see in the previous section, the fermionic free Feynman propagator in the BN 
model has a single pole which means that there are no antiparticles in the model and it coincides with the retarded 
propagator of the theory. Consequently, all closed fermion loops are zero, thus there is no self-energy correction to 
the photon propagator at zero temperature, as we already discussed it. So, we will set Q 12 = 0, therefore the closed 
fermion loops as well as the photon self-energy remain zero even at finite temperature. Another, mathematical reason, 
why we must not consider dynamical fermions - which could show up in fermion loops - is that the spin-statistics 
theorem [2] forbids a one-component dynamical fermion field. 

This means that now the exact photon propagator reads in Feynman gauge at finite temperature 


Gab,fii'{k^ — 9^vG ahij^^-iGra — ^2 


ko—^kQ-\-i€ 


And the exact photon spectral function is 


p(fc) = 27rsgn(A:o)(5(fc^). (16) 

All other propagators can be expressed using the identities in (13) and in (15). 

Now we can consider the system of self-consistent equations at finite temperatures. The derivation of the DS 
equations at non-zero temperature can be found in [7] where the CTP (Closed Time Path) formalism was applied 
(see [13]). Thus, one can express the equation for the fermion self-energy with the two-component notation as it can 
be seen in Fig. 1. In terms of analytic formulas it reads: 

'^ab{x,y) = iaae^u^ ^ / d^wd^zQac{x,w)G^''^{x,z)Ty.dcb{z]w,y), (17) 

where Ofa = (—1)“"'’^. In Fourier space it is: 

2 n 74 7 

Sa6(p) = V I Qac{p - k)G‘"^l{k)T^.acb{k]p - k,p). (18) 

This equation is the non-zero temperature equivalent of (7), the only difference here is that, since we evaluate each 
operator on a given time contour, we need to indicate them, thus we use the lower indices for this purpose. 

The vertex function in this case can be derived in a similar manner, too. The fact that the vertex function is 
proportional to is crucial again, in order to be able to apply the Ward-identities in a way we did already in the zero 
temperature computation. The derivation of Ward-identities at finite temperature is a straightforward generalisation 
of the formula that we have at zero temperature. It is easy to rewrite it in the two-component formalism, taking into 
account that to satisfy the delta functions requirement the time arguments must be on the same contour. One finds 
in Fourier space 


knKbcik-.P,q) = [^abGdi^) - ^acGdip)] {2i^f5{k+p-q). (19) 

In the BN model, because of the special property of the vertex function, it is completely determined by the fermion 
propagator in the form: 


^abc{k;p,q) 


[^abGbc (q) ^acGbc (p)] 


p—q—k 


( 20 ) 


therefore the DS equations for the fermion propagator become closed. As the third component of the system of 
equations we will use (8) again. From (8), (18) and (20) the detailed calculation of the full solution can be found 
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(a) (b) 



(c) 


FIG. 3: The fermionic spectral functions in the BN model at finite temperature. In (a) and (b) we see the 
broadening of the quasiparticle peak as the temperature and the coupling constant increased at fixed a = 0.5 and 
temperature /? = 1/T = 1, respectively. In (c) the velocity dependence of the spectral function is shown for a = 0.5. 
The r] values are rapidities, |v| = tanhry. The peak region becomes sharper with increasing rj (|v|). 


in [6]. At the end of the procedure we will get the following expression for the spectral function (defined as the 
discontinuity of the retarded propagator) in the special reference frame where u = (1, 0, 0,0): 


p{w) 


A^^/3sinQe^“/^ 
cosh{/3w) — cos a 


1 


r 




( 21 ) 


where is an a dependent normalization factor, /? is the inverse temperature and w = po — m is the variable 
(the energy shifted by the mass parameter). It can be shown that this result is completely consistent with the zero 
temperature one when one takes the T —> 0 limit [7]. In Fig. 3 we show the spectral function for different parameters 
of the temperature and the coupling constant a. By increasing the temperature the broadening of the quasiparticle 
peak can be observed and the same effect can be achieved by increasing the coupling of the theory. The halfwidth of 
the curve can be related to the inverse of the lifetime of the quasiparticle; hence increasing the temperature lowers the 
lifetime of the quasiparticle which is physically sensible. The same can be said about varying the coupling constant. 
Since the Lorentz invariance does not hold at finite temperature we need to calculate the spectral function numerically 
when one goes from the rest frame of the fermion to a different one, for details see [7]. By increasing the four-velocity 
of the fermion the halfwidth shrinks, indicating more stable quasiparticles. The finite temperature fermion propagator 
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can be retained using the spectral function in (21): 


Q{p) = J 
0 


duj p{u}) 

2tt p — u) + ie 


( 22 ) 


Since the T —>■ 0 limit of (21) gives back the spectral function that can be obtained for zero temperature (c.f. [6]), 
this limit of the finite temperature propagator also gives back the T = 0 formula (11) for this limit. 


IV. APPLICATIONS 


In this section we will present some of the ideas of the possible applications. As we have demonstrated that the 
soft photon contribution to the propagator of the electron is important both for T = 0 and T > 0 [(11) and (21)] we 
can use these expressions when we calculate processes in external fields. In such cases we can use the usual formula 
for the propagation in external field (see Fig. 4), except we use the propagator obtained from the BN model instead 
of the free one (both for T = 0 and T > 0): 


QA{xf,Xi) = Q{xf,Xi) + 


d‘^xig{xf,xi)e4-{xi)G{xi,Xi) 


+ 



Xid‘^x2g{xf,xi)e4{xi)gix2,xi)e4-ix2)g{x2,xi) + ... 


(23) 


Here we used Qa for the full propagator in the external field and Q for the BN fermion propagator. Since we have 
a simpler form of the BN propagator in momentum space we should use for our calculations: 

GAiPf,Pi) = g{Pf)i‘2TT)‘^s*{pf - Pi) + J ^^giPf)e4iPi)giPi) 

+ j j ^y^^'(P/)e4(pi)ei(K+P2)e4(p2)ei(Pi)(27r)'‘5'‘(p/-Pi-P2-Pi) + - (24) 

In this expression A'^ can be considered as an arbitrary external field. For example we can use the Coulomb potential 
generated by a static point charge —Ze: 

Ao(x) = -^, A(a:) = 0. (25) 

Or in momentum space Aq = —Ze47r/|q|. Alternatively, we can also use the screened Coulomb field: A^ = 
—Ze5^^/\x.\ exp(—jxj/?), where I is the screening length. 


+ 


+ 


+... 




FIG. 4: The BN propagator in the presence of an external field. The double line denotes the BN propagator, the 
letter ”X” is for the classical source producing the external field which are the wavy lines. The subscript A on the 
left hand side denotes that the propagator is in the external field. 


The interaction of the dressed electron with external field could have an application in the computation of processes 
like the multiphoton Compton scattering and the second order processes like the Bremsstrahlung, which are getting 
more important as the laser intensities are increasing. The main idea is similar to those of used in [14, 15] but instead 
of the Volkov propagator of the fermion the Bloch-Nordsieck propagator could be used. 


V. SUMMARY 

In this paper we discussed the infrared catastrophe in the framework of QED and its resolution using the Bloch- 
Nordsieck summation of radiative corrections. We sketched an alternative derivation of the Bloch-Nordsieck propa¬ 
gator with the aid of the Ward-identities; the full derivation is presented in [6] and for finite temperatures in [7]. The 



Bloch-Nordsieck propagator describes the propagation of the dressed electron in the deep infrared regime of the QED, 
where the soft photon contributions can be completely summed up. Speculation of the possible applications of the 
results are given in the last section. The detailed calculations of these processes are topic of current research projects. 
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